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I. INTRODUCTION 



Standard Random Phase Approximation (s-RPA) is 
one of the most popular many body approches known. It 
was invented in condensed matter physics (see e.g. 0) 
and has subsequently spread to almost all branches of 
hysics, including atomic physics 0, molecular physics 
, plasma physics Q , relativistic field theory |5l| , nuclear 
physics 0, and many more. The definition of s-RPA is 
not uniform, depending on whether exchange is included 
or not. We understand it, e.g. as in nuclear physics Q, 
as the small amplitude limit of time dependent Hartree- 
Fock theory (TDHF) and therfore with exchange. Rs 
popularity probably stems from its conceptual simplicity, 
its numerical tractability (in spite of some serious prob- 
lems in finite size systems), and most of all its well be- 
haved properties conserving full fiUment of conservation 
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laws (Ward identies), Goldstone theorem, and restoration 
of spontaneously broken symmetries. Though there exist 
respectable general theories (see e.g. 0, 0|), any prac- 
tical attempt to go beyond this basic HF-RPA scheme 
conserving these properties turned out to be technically 
extremely demanding and no well accepted general and 
practical extension has emerged so far. Nevertheless, 
standard RPA has also quite serious shortcomings and 
it is desirable to overcome them. One of the most promi- 
nent is its violation of the Pauli principle, often para- 
phrased as the "quasiboson approximation". It is most 
critical for only moderately collective modes or when the 
self interaction of the gas of quantum fluctuations be- 
comes important as in ultra small finite quantum sys- 
tems. Since a couple of years two of the present authors 
and collaborators have been working on a non linear ex- 
tension of RPA Q which has shown surprisingly accurate 
results in a number of non-trivial models jlOi | . It is called 
Self -Consistent -RPA (SCRPA) and can be obtained 
from minimising an energy weighted sum rule. Therefore 
s- RPA which is perturbative in the sense that it sums 
a certain class of diagrams (the bubbles) is upgraded in 
SCRPA to a nonpcrturbative variational theory though it 
is in general not of the Raleigh-Ritz type. A strong bonus 
of this extension of s-RPA is that it generally preserves 
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its positive features as conservation laws and restoration 
of symmetries as well as numerical tractability, since it 
leads to equations of the Schrodinger type jlJl- In this 
paper we want to apply this theory to the Hubbard model 
for strongly correlated electrons. Because of its necessar- 
ily increased numerical complexity over s-RPA, we first 
want to consider finite clusters in reduced dimensions. 
Before going into the details, let us very briefly repeat 
the main ideas of SCRPA. 

One way of presentation is to outline its strong anal- 
ogy with the Hartree-Fock-Bgoliubov (HFB) approach to 
interacting boson fields b\ b. The HFB canonical trans- 
formation reads 



with 



Vi,i, hi 



(1) 



The amplitudes m, v can be determined [12| from min- 
imising the following mean energy (energy weighted sum 
rule) 



(0|[(7.,[i?,4]]|0) 



01 



ql |0 



(2) 



where H is the usual many body Hamiltonian with two 
body interactions and the groundstate |0) is supposed to 
be the vacuum to the quasiboson operators g^, i.e. 



g.|0) = 



(3) 



With this scheme and the usual orthonormalisation con- 
ditions for the amplitudes u, v, which allows the inver- 
sion of 1^, one derives standard HFB theory .6] with no 
need to construct |0) explicitly. Of course, in this way 
the fact that HFB is a Raleigh-Ritz variational theory 
is not manifest but the scheme has the advantage to be 
physically transparent and to lead to the final equations 
with a minimum of mathematical effort. 

For SCRPA we follow exactly the same route. We re- 
place in Q the ideal boson operators by fermion pair op- 
erators of the particle-hole (p/i) type and form an ansatz 
for a general transformation of ph -Fermion pairs 



ph 



(4) 



with — QI\0) an excited state of the spectrum. In 
analogy with we minimise a mean excitation energy 



(0| [Q,, [H, Ql]] |0) 



0| Q., Ql |0 



(5) 



with 1 0) , in analogy with Q , the vacuum to the operators 
Q„, i.e. 

Qu\0) - (6) 
and arrive at equations of the usual RPA type 



-B* -A* j ^ " [ y 



(7) 



A 



ph,p' h' 



al ap , 



H, aj,, ah' 



- npy/Hh' - Tip' 



ph.p'h' 



al Qp , 



H, a\, ap' 



(8) 



Here we supposed to work in a single particle basis which 
diagonalises the density matrix (natural orbits) : 



\al ak'\0) = rikSkk' 



(9) 



and therefore the n^'s are the occupation numbers. For 
H with a two body interaction, ((SJ only contains corre- 
lation functions of the (a^ a) and {a^ a a^ a) types and, 
since Q admits the usual RPA orthonormalisation rela- 
tions for the amplitudes X, y 6], the relation Q can be 
inverted and with ((HJ the correlation functions in (jSJ be 
expressed hy X, y. 

However, to be complete, occupation numbers Uk = 
(0|a|. afe|0) and two body correlation functions with other 
index combinations than two times particle and two times 
hole need extra considerations. That will be done in the 
main text. This is, in short, the SCRPA scheme which, 
as HFB theory, is obviously non linear, since the elements 
A and B in ^ become functionals of the X and y am- 
plitudes. We want to point out that no bosonisation of 
Fermion pairs is operated at any stage of the theory. 

We want to apply this scheme to the Hubbard model 
of strongly correlated electrons which is one of the most 
wide spread models to investigate strong electron corre- 
lations and high Tc superconductivity. Its Hamiltonian 
is given by 



H 



E 

<ij><y 



^E^ 



(10) 



tion operators at site 'i' and the hi^ = Cia are the 
number operators for electrons at site 'i' with spin pro- 
jection a. As usual t is the nearest neighbour hopping 
integral and U the on site coulomb matrix element. In 
this exploratory work, we will limit ourselves to the sim- 
plest cases possible, i.e. we will consider closed chains 
in one dimension with increasing number of sites at half 
filling, starting with the 2 -sites problem. It will turn 
out that the next case of 4 -sites is a configuration with 
degeneracies which cause problems in SCRPA, as do all 
4n [n = 1,2,3,...) configurations in ID. We therfore 
will postpone the treatment of these cases to future work 
and directly jump to the case of 6 -sites and only shortly 
outline at the end why the 4 -sites case is unfavorable 
and how the problem can eventually be cured. In this 
work we will stop with the 6 -sites case considering it as 
sufficiently general to be able to extrapolate to the more 
electron case. In this way one may hope to approach the 
thermodynamic limit in increasing the number of sites 
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as much as possible. Let us mention that an earher at- 
tempt to solve SCRPA in ID in the thermodynamic limit 
in a strongly simplified version of SCRPA, the so-called 
renormalised RPA (r-RPA), produced interesting results 

m 

In detail our paper is organised as follows : in section 
(|IT| we present the two sites case with its exact solution. 
In section (|III|I we outline the 6 -sites case with a detailed 
discussion of the results and in section (jIV|) we present 
the difficulties encountered in the 4 -sites case and how, 
eventually, one can overcome them. Finally in section (|V|) 
we give our conclusions together with some perspectives 
of this work. 



II. 



THE TWO SITES PROBLEM 



In this section we will apply the general formalism of 
SCRPA outlined in the introduction to the two sites prob- 
lem at half filling, i.e. two electrons with periodic bound- 
ary conditions. This case may seem trivial, the fact, how- 
ever, is that such popular many body approximations as 
standard RPA (s-RPA), GW [lil, Gutzwiller wave func- 
tion [13, the TPSC approch by Vilk and Tremblay [Tfil |. 
etc.. do not yield very convincing results in this study 
case, whereas it has recently been shown that SCRPA 
solves two body problems exactly 0, 0, 0| . We again 
will briefiy demonstrate this here for the two sites prob- 
lem. 

First we will transform 11(J|) into momentum space. 
With the usual transformation to plane waves 



■j this leads to the standard expression 



for a zero range two body interaction 



H = 



k,cr 

u 



k,p,q,a" 



where fi-^^ — o,^^ is the occupation number operator 
of the mode (k, a) and the single particle energies are 

D 

given by = ~2t ^ cos (kd) with the lattice spacing 



set to unity. 

For our further considerations it is convenient to trans- 
form to HF quasiparticles operators via (we switch 
to II?) 



a/i,, 



(12) 



where h and p are momenta below and above the Fermi 
momentum, respectively, so that bk,cr\HF) = for all 
k where \HF) is the Hartree-Fock groundstate in the 
plane wave basis. For the two sites problem with pe- 
riodic boundary conditions we then write, after normal 



ordering, the Hamiltonian in the following way 

H = Hhf + Hq=0 + Hq^T^ (13) 

with 



Hhf — E 



HF 



El 



El 



u 



U 



Hq=Q = ^(^fe2,T-"fci,T)("fe2a-^fci,i) (15) 



|(j-+J+)(j; + J+) 



(16) 



and J„ = &i,o-62,<T, J+ = [Ja)^ 1 "-fci,<T = b\^^hi^a 
where we introduced the abreviation "1" and "2" for the 
two momenta fci = and fc2 = — tt of the system, re- 
spetively. The HF groundstate is \HF) = bi^^bi^i\vac) 
and the corresponding energy is given by 

E^^ = {HF\H\HF) = -2 1 + ^ (17) 

The RPA excitation operator corresponding to can, 
because of rotational invariance in spin space, be sepa- 
reted according to spin singlet (S = 0, charge) and spin 
triplet {S = 1) excitations. The latter still can be di- 
vided into spin longitudinal (S* = 1, = 0) and spin 
transverse {S = 1, rris = ±1) excitations. Let us first 
consider the charge and spin longitudinal sectors. For 
later convenience we will not separate them and write 
for the corresponding RPA operator 



Ql = A-f K+ + XI K+ - y\ K- - y'l K 



(18) 



where ^ J^/ ^1 - (M^), M„ = ni„ + n2a , and 
the mean values (. . .) are always taken with respect to 
the RPA vacuum 

Q,\RPA) = . (19) 
Because of the orthonormality relations 



[xix'^,.-yiyi,)^8. 



/ J \ (J a 



J2 iKy:'-y:K') = o, 



(20) 



one can invert (|18|l to obtain 

J- = v/1 - (M^) J2 ( -^^ Q-^+yaQt) , 

Ji = • (21) 
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The operators and 1 — M„ form a SU{2) algebra of 
spin -i operators and, therfore, using the casimir relation 
we obtain 



2J+J- 



(22) 



In this way we can calculate with H19|) the following ex- 
pectation values 



{J-> Jt ) = v/(i-M.-)(i^M.) J2 K, 

V 

{Ji> Ji ) = ^{i-M,.){i-M,) J2 y^' K, 

V 

Ja ) = y/{l-M„,){l-M,) J2 K' (23) 



with 



{Ma 



2 E 13^.1' 

1 + 2 E |3;-| 



(24) 



We will see that in order to close the system of SCRPA 
equations, expectation values {McrM^') will also be 
needed. It is easy to see that we have 



and 



MaMa' - 4 44, Ja' Ja {(T ^ o') 



(25) 



(26) 



With (|21() the expectation value of H26() gives 
{MaMa') = 4(1 - (A/^))(l - {Ma')) 

■ E E y^yay^^y:^ {q.q.^q^qI) ■ 



(27) 



For the calculation of the correlation functions which ap- 
pear on the right hand side of H27() one commutes the de- 
structors Qi, to the right and uses © , yielding again cor- 
relation functions {M^Ma')- One then obtains a closed 
linear system of equations for the latters. Details are 
given in Appendix 

The SCRPA matrix elements can be expressed in the 
following way 



= 2f + S|,| 

= 2t + Si,i 

= 

- Bu (28) 





K-, 


H, K+ 


Ad = ( 






-^T.i = ( 






^i.T = ( 




H, K+ 



B 



T,T 



K 



= u. 



Bu = 



' l-(Mi) 
1 - (Mt) 



E i^iyi+'^iXi) 



'1- (Aft) 



(Mi) 



E i^ry: 



y^yn 



U ((l-AfT)(l-Mi)) 
2 ^(l-(AfT))(l-(M,)) 



B 



H, K 



B 



(29) 



With our previous relations (|23|l . (|24|l and (|27|l we can 
entirely express the elements of H28(l and (|29|l by the RPA 
-amplitudes and therefore we have a completely closed 
system of equation for the amplitudes A", y. With the 
orthonormality relations (|20|l we further more have 



Aa 



Ad 
B^d- 



--Ad 

Bid- 



B' . 



(30) 



and, therefore, the SCRPA equation can be written in 
the following form 



/ A A' B B' 

A' A B' B 

-B -B' -A -N 

\-B' -B -A' -A 



The system 



has 



y^ 
\yi J 

the 



yyi J 



two 



(31) 
positive 

and 



roots £i = \J{A- A') - {B - B') 

82 = ^JiA + A')"^ - [B + B'f. The SCRPA equa- 
tion (|31|) can be solved numerically by iteration leading, 
as expected, to the exact result. This latter fact can also 
be seen analytically in noticing that by symmetry 



yl = -yl = y,. 



, = yl = y-h ■ (32) 



Therefore the 4x4 equation 131|) decouples into two 2x2 
equations corresponding to charge {ch) and spin {sp). 
Then we see that the exact groundstate wave function 
which contains only up to 2p — 2h excitations 



|0) cx (l + rf J+ J+) \HF 



(33) 



is the exact vacuum to the RPA operators, i.e 
Qch{sp)\RPA = under the condition that 



= tan {(j)) . 



(34) 



ch{sp) 



We therfore can express the SCRPA equations by the 
single parameter (j) and obtain the solution analytically 
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(up to the solution a nonlinear equation for (p). The 
solution agrees for all quantities with the exact result. 
For example the groundstate energy is given by 

eSCRPA ^ _2 <cos(2 fli) + ^ (1 - s«n(2 4>)) ■ (35) 

This expression can either be derived directly from (H) 
using and or one uses a generalisation of the 
standard RPA expression for the groundstate energy : 



jT^SCRPA 



(36) 

It is straightforward to verify that expressions (|35|l and 
are identical. 

The standard RPA expression are recovered from H31|) 
in replacing in all expectation values the RPA ground- 
state by the uncorrelated HF determinant. In Fig^we 
compare the standard RPA with SCRPA and exact re- 
sults for the excitation energies and in Fig|21 the cor- 
responding groundstate energies together with the HF- 
values are shown. From these figures one should espe- 
cially appreciate the long way SCRPA has gone from s- 
RPA to recover the exact result. For instance it is clearly 
seen that the instability of s-RPA atU — 2 is, as expected 
for such a small system, an artefact and is completely 
washed out by the self consistent treatment of quantum 
fluctuation contained in the SCRPA approach. 



ph -RPA Standard 

H ^ ph -SCRPA 

Exact 




FIG. 1: Excitation energies of the standard RPA (dashed 
lines), SCRPA (crosses) and exact solution (solid lines) as a 
function of U in the channels of charge (ch) and longitudinal 
spin (sp) for the 2 -sites case. 



Without explicit demonstration let us also mention 
that SCRPA in the spin transverse channel with Ql, = 

as well as in the particle-particle channel with = 
X ~ y ^li ^iT ^^^o gives the exact solution for the 

two sites problem. How the pp-SCRPA works can be seen 



m ref[l3 where for the pairing problem the two particle 
problem is also solved exactly. 



-0.5 



E„Jt 



GS" -1 




HF 

ph-SCRPA 

— Exact 

— ph -RPA standard 



4 
U/t 



FIG. 2: Groundstate energy in HF (dot-dashed line), stan- 
dard RPA (dashed line) , SCRPA (crosses) and exact solution 
(solid line) as a function of U in the charge and longitudinal 
spin responses for the 2 -sites case. 

The fact that SCRPA solves the two sites problem ex- 
actly is non trivial, since other well known many body ap- 
proaches [Til llSl llq. as already mentioned, so far failed 
to obtain this limit correctly. 



III. 



THE SIX SITES PROBLEM 



After this positive experience with the two sites prob- 
lem we next will consider the 1-dimensional 6 -sites case, 
as for the 4 -sites case problems appear needing particu- 
lar considerations to be outlined in section llVl We again 
consider the plane wave transformation explained in sec- 
tion ^ with the corresponding Hamiltonian in momen- 
tum space In the first Brillouin zone — tt < k < tt 
we have for — 6 the following wave numbers 



fci 

ki — 



TT 

3 



27r 



, k2 

T ' 3 

With the HF transformation 



k. = — 



(37) 



(38) 



such that bk^a\HF) — for all fc, we can write the hamil- 
tonian in the following way (normal order with respect 
to b\ b) 



H — Hhf 
where 



191=0 



l9l = f 



H\q\=Tr (39) 



6 



H 



HF 



H\q\=0 



3 



J. 



5T,3T 



S. 



6i,4i 



5", 



5i,6i 



- s 



s 



G 



5T,4T 



'^3T,2T 



+ 



U,2i 



J- 



3i,5i 



6T,2T 



5- 



4i,5i 



s. 



G 



Ji 



T,6T 



2i,3i 



51M 



Jl 



lAl 



J. 



6i,3i 



J. 



21.61 



'3T.4T 



'2i,5i 



'3i.4i 



CC 



(40a) 
(40b) 

(40c) 

(40d) 
(40e) 



with the following definition of operators 



Jph.a - bh "b 



Jph.a - (Jph.a) (41) 



and 



£l 
£4 
G 



87+ ^t/, 

-27+|, 
7/ 

£5=t+2 > 

U 

"6" 



£2 = £3 = 
£6 = 27 



7/ 



7/ 
2"' 



(42) 



The level scheme is shown in FigEl The hole states are 
k=k,=-7r A I 



V 



A 



A 



k=k^=27c/3 





/i 




/i 




i 




i 







/i 








i 





FIG. 3: Excitation spectrum of HF at [/ = for the chain 
with 6 -sites at half filling and projection of spin nis = 0. The 
occupied states are represented by the full arrows and those 
not occupied are represented by the dashed arrows. 



\i\ = ¥ 


\q\ = TT 


l9l = f 


51 ^ 951 = 


61 qui ~ —71" 


42 ^ 942 = -h^ 


41 ^ 941 = +^ 


52 952 = — TT 


53 ^ 953 = 


62 ^ = +^ 


43 943 = -fTT 




63^963 = -^ 







TABLE I: The various momentum transfers in the 6 -sites 
case. 



(43) 



The HF groundstate is 

\HF) = ai_^a|^^4^^4_^4^^4 J-) . 

We see that the Hamiltonian for 6 -sites has largely 
the same structure as the one for 2 -sites. It is only 
augmented by H\q\^:^ + i?|q|^2x which contains the 5- 
operators on which we will comment below. 

There are three different absolute values of momen- 
tum transfers as shown in Table Since the momentum 
transfer |g| is a good quantum number, the RPA equa- 
tions are block diagonal and can be written down for each 
\q\ -value separately. For exemple for |g| = ^ we have the 
following RPA operator for charge and longitudinal spin 
excitations 



'^21, 4T -'^41,21 



'^3't,5T -^5T:3T + '^3'j.,5i -^51, 3i 



3^2T,4T -^2T,4T 
^3T,5T -^51,31 



3^2i,4i -^2i,4i 
^3J„5i-^3i,5i 



where 



and 



K: 



p(T,hcr 



pa, ha 



(44) 



(45) 



labeled h = {1,2,3} and the particle states p = {4, 5, 6}. 



M, 



pa.ha 



(46) 
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We write this RPA operator in short hand notation as 
again with the properties 



{xr Jt - y- jn (47) 



Q.|0) = . 



(48a) 
(48b) 



The matrix elements in the SCRPA equation 

-B* -A* j \ y j ^ [ yi^ ^ are then of 
the form 



Aii' — 



Bu' 



V(1-(M,,))(1-(M,)) ' 
V(1-(M.'))(1-(A/.)) 



(49a) 
(49b) 



Since the SCRPA equations have the same mathemat- 
ical structure as standard RPA, one also has equivalent 

orthonormahty relations J2 (K^K"' ~ yi y^') = ^I'u- > 

etc.. , in analogy to eqs. H20() of the 2 -sites case. This al- 
lows to invert H47|l and to calculate the expectation values 
which will appear in (|49a|l and Il49b|) in complete analogy 

to HSli. 

The missing expectation values {Mi) can be expressed 
by the X,y amphtudes in observing that jf, and Jf — 
i(Mi — 1) form, as in the 2 -sites case, an SU2 Lie algebra 
for spin-i particles. Using the Casimir relation one again 
obtains A'L = 2 Jt J7 and thus 



2 



i/|2 



1 + 2 Ei3^ri 



(50) 



We also will need expectation values of 



M,M, =4J+J-J+J^ 



for i 7^ j , 



(for MiMi — 2 Mi we can use (l.'iOll '). Those can again 
be calculated following the same procedure as outlined 
in 1211) and Appendix 



In order to solve the SCRPA equations we now practi- 
cally have prepared all we need. Nontheless, at this point 
we have to discuss a limitation of our RPA ansatz 144|) 
which is not absolutely necessary but which turned out 
to be convenient for numerical reasons. The fact is that 
our RPA ansatz is restricted to ph and hp configurations, 
as this is also the case in standard RPA. In the latter case 
this is a strict consequence of the use of HP occupation 
numbers rip, n'^ with values zero or one, respectively. In 
the SCRPA case with a correlated groundstate the oc- 
cupation numbers are different from zero and one and 
a priori there is no formal reason not to include into 
the RPA operator also pp and hh configuration of the 
form aj,ap' = bj,bp' and a\ah' = —b]^,bh- Such terms are 
usually called scattering or anomalous terms T^. With 
rounded occupation numbers the SCRPA equations (at 
T = 0) are formally and mathematically equivalent to 
standard RPA equations at finite temperature where also 
pp and hh components are to be included, in prin ciple 
|18| . The inclusion of those scattering terms 

(the 

S-terms in ()39|l ') usually is of little quantitative conse- 
quence entails, however, the important formal prop- 
erty that, as for standard RPA, the energy weighted sum 
rule is fuUfiUed exactly 0, Inspite of this desir- 

able feature, we had to refrain from the inclusion of the 
scattering configurations in this work because the factors 
y/l — {Mi) by which the SCRPA matrix is divided (see 
eqs. H49a|) and (|49b|) '). can become very small in these 
cases and this perturbed the convergence process of the 
iterative solution of the SCRPA equations. Though we 
do not exclude that a more adequate numerical proce- 
dure could be found to stabilise the iteration cycle, we 
decided to postpone such an investigation, because, as 
already mentioned and as will be shown later, the influ- 
ence of the scattering terms is, as found already in other 
studies very small. We will shortly come back to 
this discussion when presenting the results for the en- 
ergy weighted sum rule below. As a consequence and 
for consistency we then also will have to disregard the 
5-terms of the Hamiltonian (remember that also in stan- 
dard RPA these terms do not contribute). Under these 
conditions we then obtain a completely closed system of 
SCRPA equations. For completeness we give some ex- 
amples of SCRPA matrix elements which correspond to 
the ansatz for 1^1 = ^ 
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A2A 



(1 - (Af24,T» 



(A- 



(1 - (A/24,T»" 



'2i,4i 



4T,2T 



V(l - (Af24,i» (1 - (M24,T» 

g|((1 - Af24,T) (1 ~ Af24,i)) + ((4.1T - ^6^21) (^il 

'•^3^.41 - ^2^.5|)>| {(1 - (^'^24a» (1 - (M24,T»r 



U,4i "^21,61 



+ ^("^4T,3T "^51,21; V"^3i,4i ^2i,5i 



(51a) 



(51b) 



The other matrix elements can be elaborated along the 
same lines. Of course in the approximation where the 
expectation values in (|51a|l and Ij51bp are evaluated with 
the HF groundstate the usual matrix elements of stan- 
dard RPA are recovered. We should also mention that in 
expressions (|51all and 151bp expectation values as for ex- 
ample ( Jj^ 4| J^^ which involve momentum transfers 
other than the one under consideration (jgaj = ^ in the 
specific example) must be discarded. That this implicit 
channel coupling cannot be taken into account without 
deteriorating the quality of the SCRPA solutions is an 
empirical law which has been established quite sometime 
ago • It is part of the decoupling scheme and it is in- 
tuitively understandable that, since each channel is sum- 
ming specific correlations, one can not mix the channels 
implicitly without perturbing the balance of the minimi- 
sation procedure which is done channel by channel. It 
can also be noticed that, neglecting the S'-terms in H, 
the channel coupling disappears. 



We here give for the transfer |g| = ^ the totality of the 
elements of matrix SCRPA A and B just as it is was used 
in the numerical calculation. For others transfers there 
will be analogous expressions. Indeed with the following 
abbreviations 



1 = (2 T,4 T) 
3 = (3 T,5 T) 



2=(2i,4i) 
4=(3i,5i) 



the elements of matrices A and B are given by 

i _ ^ ^ ^^ (-^21,41 (-^31,51 + '^44,24)) 

A.i - e4-e2-2G YZJ^^^) ' 

^ C ((1-M24,T)(1-M24a)) 

V(l-(Af24,i))(l-(M24,T)) ' 

-^3,1 ~ -^4,1 — Az,2 ~ Aa^ = , 
^ ^ ^ r, ^ ^(^3T,5T + '^4T,2t)'^2X,4x) 

-42,2 - e4-62-2G Y^^) ' 

A _ ^ ^ (-^31,51 (-^2J„4i + ■^5'l,3J,)) 

y^3,3 - £5 - ea - G — , 

1 - {Ahb,] ) 

^ ^ C ((l-Af35.T)(l-A-^35a)) 

v/(l-(A/35a))(l-(A/35.T» ' 



^4,4 — £5 ~ £3 

K _ ^^ ('^2T,4T('^2i,4J, + '^5J„3i)) 

1 - (M24.T) ' 

^2,1 = ^3,1 = 04,2 = 04,3 = , 

g ^ ^ ((l-Af24,T)(l-Af35,i)) 

v/(l-(Af35,i))(l-(A^24,T)) ' 

H _ „^(('^2T.4T + "^5T,3T)"^2i,4i) 

^^'^ ^ '''' T^W^) ' 

g ^ ^ ((l-Af35.T)(l-Af24,i)) 

V(l-(A/24,i))(l-(Af35,T» ' 

R _ r,^^'^3T,5T("^3i,5i + •^4J„2i)) 

- 1 - (M35,t) ' 



04 4 = -2G 



(("^31,51 + "^4T,2T)"^3i,5i) 



(52b) 
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Let us add that the matrices A and B are symmetric 
and that the expectation values (. . .) in H52a|) and Ij52b|l 
can be expressed in an analogous way as the expectation 
values (|23|l and (|27|l by the amplitudes X, y. 

The structure of the self consistent matrix elements 
(|52a|l and (j52bp is also quite transparent : the bare in- 
teraction which survives in the limit of standard RPA is 
renormalised, i.e. screened, by two body correlation func- 
tions which are calculated self consistently. The general 
structure of the scheme is in a way similar to the one 
proposed by Tremblay and coworkers p^ . however, the 
details of the expressions and the spirit of derivation are 
different. One can also interpret our theory as a mean 
field theory of quantum fluctuations as this was done in 
refs.0 



lql=7i: 



lql=27c/3 




0.5 




pli -RPA standard 

Exact 

^ ph -SCRPA 



0.5 1 1.5 2 2.5 3 3.; 

u/t 



FIG. 4: Energies of excited states in standard RPA, SCRPA, 
and exact cases as a function of U for 6 -sites with spin pro- 
jection rris — and for |q| — n. States of the charge response 
and those of the longitudinal spin response are denoted by ch 
and ap, respectively. 

Let US now come to the presentation of the results. In 
Fig. 10] [S] and O we display the excitation energies in 
the three channels \q\ = tt, ^ and ^ as a fonction of 
U/t. The exact values are given by the continuous lines, 
the SCRPA ones by crosses and the ones corresponding 
to standard RPA by the broken lines. We see that in 
all three cases SCRPA results are excellent and strongly 
improve over standard RPA. As expected, this is par- 
ticularly important at the phase transition points where 
the lowest root of standard RPA goes to zero, indicating 
the onset of a staggered magnetisation on the mean field 
level. It is particularly interesting that SCRPA allows 
to go beyond the mean field instability point. However, 
contrary to the two sites case where SCRPA, in the plane 
wave basis, solved the model for all values of U ^ here at 
some values U > Ucr the system "feels" the phase transi- 
tion and SCRPA stops to converge and also deteriorates 
in quality. Up to these values of U SCRPA shows very 



FIG. 5; Same as Fig. Hbut for \q\^'^. 



good agreement with the exact solution and in particular 
it completely smears the sharp phase transition point of 
standard RPA which is an artefact of the linearisation. 

lql=Jl/3 



■ ph -RPA standard 

■ Exact 

^ ph -SCRPA 




FIG. 6: Same as Fig. Hbut for |ql = f 
In Fig. [7|we show the groundstate energy (see eq. (|36|l ') 



rpSCRPA 

^0 



(53) 



as a function of U. In addition to exact, SCRPA, and 
s-RPA values we also show the HF energy. Again we see 
that SCRPA is in excellent agreement with the exact so- 
lution. Standard RPA is also good for low values of U but 
strongly deteriorates close to the lowest phase transition 
point which occurs in the \q\ = tt channel at U — 
The HF energies, on the contrary, deviate quite strongly 
from the exact values. 
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ph -SCRPA 
HF 




0.006 



0.004 



0.002 



(p,h)=(2,4) 

(p,h)=(3,5) 




FIG. 7: Energy of groundstate in HF, standard RPA, SCRPA 
and exact cases as a function of U for 6 -sites with spin pro- 
jection TUs — 0. 



FIG. 8: The ratio, r (eq. I54|l . as function of the interaction U 
for the ph excitations (2,4) and (3,5) in the channel |g| = ^. 



with 

The reader certainly has remarked that our RPA 
ansatz (|44|) has so far not separated charge and spin ex- 
citations. In the 2 -sites problem this was automatically ^ _ _ ^^-j j^j/ o)|^ 
and exactly the case. However, here, since we did not 

consider the S* -operators neither in the Hamiltonian nor x^/ w/ \/^ \ 

in the RPA operator, spin symmetry is violated. On the = " ^o) mQ\q\,^ F\0}\ (56a) 

other hand this permits to evaluated the importance of "'■'''I 

the S -operators. Normally the eigenvectors of the RPA = ^ _ j^^'j |^0| [Q\q\ v, F\ |0)p 

matrix should be such that for charge (eft.) excitations ^ 

the operators Jth^+Jth\ ^-^d J~ht+J'h, can be fac- _ _ , 

tored whereas for spin (sp) excitations the combinations 2-^ ^ ' v ^ ^"i ~yi ) \ 

Jph^ - Jphi and J-^^ - J-^^ hold. Because of our vi- -■\<i\ 
olation of spin symmetry this factorisation is not exact. ^ _ iQ\\F \H F]\\(S) 

To have a measure of this violation we plot in Fig. |H1 the 2 ' 

ratio 



(54) 



^ ^J\ - - M,, (A,,' - S,,,0(56b) 



I'^ph^l + I'^phil with 



For exact spin symmetry r should be zero. From Fig. |S1 
we see that the violation is on the level of a fraction of 
one percent. This, therefore justifies, a posteriori having 
neglected the scattering terms {S -terms) in the Hamil- 
tonian and RPA operator. A further indication that S - 
terms are not important comes from the energy weighted 
sum rule. We know that the sum rule including the S 
-terms is fullfiUcd in SCRPA 13, 19] . However, neglect- 
ing them gives a slight violation. Considering the exact 
relation 



L = R (55) 



we trace in Fig. El the ratio ^ = . Again we see 
that the violation is on the level of a fraction of one per- 
cent, confirming the very small infiuence of the scattering 
terms. 

A further quantity which crucially tests the ground 
state correlations are the occupation numbers. We have 
no direct access to them, however, we will use the so- 
called Catara approximation for their evaluation 'ST] : 
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0.006 




FIG. 9: The ratio, ^ — ^^^^ , of the energy weighted sum rule 
in the charge response for the 6-sites case. 



riper — {flpa) — 'y^.jJ^h.a '^vh.a) 
h 

= Y.^l-{Mp,,.))Y.\y;naW (58a) 

h V 
P 

= 1 - ^ (1 - {Mpna))Y. \yvhA^ (58b) 

p V 

We show these quantities in Figs. llOlandlTTlin compari- 
son with the exact values and the ones of standard RPA. 
We again see the excellent performance of SCRPA. 

Concluding this section we can say that the expecta- 
tion we had from the 2 -sites case, with its exact solution, 
have very satisfactorily also been fuUfilled in the 6 -sites 
case. However, in spite of the very good performance of 
SCRPA, there is the limitation that SCRPA, in the sym- 
metry conserving basis of plane waves used here, can not 
be employed in the strong U limit. One also may won- 
der how the extension to cases with sites number 2 + An 
with n > 1 works. For such cases it does not make sense 
any more to elaborate the Hamiltonian in its detailed 
form as given in ea. (|40|l . This explicit expression was 
only given to make clear the detailed internal structure 
of the approach for a definite example. In the general 
case with many sites one would just take the form Hll|) 
of the Hamiltonian, calculate the double commutators as 
needed in ((SJ) and then express the resulting correlation 
functions by the X- and y- amplitudes. That such a 
program is feasable in terms of analytic work and nu- 
merical execution was demonstrated in our earlier work 
on the multilevel pairing model [Tol | where cases up to 
hundred levels were treated. However, this number was 
not considered of an upper limit. Though the present 
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0.5 1 1.5 2 2.5 3 3.5 

U/t 

FIG. 10: Occupation numbers as function of the interaction 
U for various values of the momenta fc for states above the 
Fermi level. For each approximation, s-RPA and SCRPA, the 
occupations numbers are represented in increasing order like 
k (— TT, — ^). Let us notice that the modes k — ^ and 
k = are degenerate. 




U/t 

FIG. 11: Occupation numbers as function of the interaction 
U for various values of the momenta k for the holes states. 
For each approximation, s-RPA and SCRPA and exact solu- 
tion, the occupation numbers are represented like fc = 0, ^, 
— -1. Let us notice that the modes fc = f and k — are 
degenerate. 



model is slightly more complicated, we think that a gen- 
eralisation to the case of many sites is perfectly possible. 
It needs, however, some investment which is planned for 
the future. This also concerns the D — 2 case. An- 
other question to ask is whether the degradation of the 
SCRPA results going from the = 2 to the N = 6 case 
does not go on considering N — 10, 14, etc? One again 
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may cite the experience with the multilevel pairing model 
[iQi] where also the N — 2 case turned out to be exact in 
SCRPA but not the other cases. However, all iV > 2 cases 
showed more or less the same degrees of accuracy: excel- 
lent results of SCRPA up to the phase transtion point 
and deterioration beyond. Since this behavior has also 
been found in simpler models |12| we think that this is 
a generic feature of SCRPA and that this behavior will 
also translate to the case of the present model. 

Another problem for further work is how to continue 
the present theory into the strong coupling regime. Of 
course, there exists the possibility to perform SCRPA 
in the symmetry broken basis, but details and how to 
match with the summetry unbroken phase must still be 
worked out. Also the inclusion of higher order operators, 
as will shortly be discussed in the next section, may be 
an intersting direction in this respect. 

IV. FOUR SITES PROBLEM 
A. The symmetry unbroken case 

The problem of the 4 -sites case is easily located in 
regarding the level scheme of Fig. ^1 (see also ref 
deeling with the attractive Hubbard model in ID). We 

k=-jt I A 

V I 



1 A 






V ; \ 









1 











FIG. 12: Level spectrum for U — for the half filled chain 
with four sites with spin projection nia — 0. The occupied 
states are represented by the full arrows and those not occu- 
pied are represented by the dashed arrows. 

see that the Fermi energy coincides with the second 
level which is half filled. The uncorrelated groundstate 
is therefore degenerate and excitations with momentum 
transfer |g| = tt cost no energy. On the other hand for 
excitations with l^l = § there is no problem. The corre- 
sponding RPA operator is given by 

— 3^13 I ^^^13 I — 3^24, T ^24, T 

- 3^l\i^r3a- ^24,1 ^2^4,1 (59) 

In Fig. Uniwe show the results of s-RPA and SCRPA, to- 
gether with the exact solution. We see that the lower ex- 
citation is still very well reproduced by SCRPA, whereas 
for the second excited state SCRPA only reduces the dif- 
ference of s-RPA to exact by half. The real problem 



shows up for the transfer \q\ = vr. The corresponding 
operator is 

+ -^2^3,1 -^32, T '^23, i ^32, i 

- yki ^23.T " yki ^23,i (60) 

The standard RPA produces a doubly degenerate zero 
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FIG. 13: Energies of excited states with standard RPA, 
SCRPA, and exact solution for four sites with spin projec- 
tion rus = and for \q\ = ^ in the symmetry unbroken basis. 



Iql=7t 



£/t 




FIG. 14: Energies of excited states with standard RPA, 
SCRPA, and exact solution for four sites with spin projec- 
tion rris — and for \q\ = vr in the symmetry unbroken basis. 

mode independent of U as seen on Fig. 1141 As compared 
with the exact solution, we see that these two zero modes 
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approximate two very low lying exact solutions. Unfortu- 
nately, because of these modes at low energy the SCRPA 
could not be stabilised. The only possibility consisted in 
excluding the components K^^ | ^'^^ ^32 x the RPA 
operator. Then the self consistency was achieved without 
problem and the result is shown in Fig. 1141 The result 
of SCRPA is half way in between s-RPA and the exact 
solution. On the other hand, because of the omission of 
the two lower states, the groundstate energy can not cor- 
rectly be calculated in SCRPA. Therefore, for the 4 -sites 
problem in the symmetry unbroken basis (plane waves), 
the SCRPA cannot fully account for the situation. 



B. Symmetry broken basis 

An analysis of the HP solution shows that, as soon as 
U ^ 0, the plane wave state becomes unstable and the 
system prefers a staggered magnetisation. The general 
HP transformation can be written as 



0.75 
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FIG. 15: Occupation numbers for site 1, ni,| et ni,|, as a 
function of interaction U in the symmetry broken basis. 
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±i. We also have two other excitation opera- 



and 



K+ 

T,4T -^414! 



^11 



3^2T,3T -^2T,3T 



lAl ^4i,U 



" '^21,31 -^31, 2i 
^2i,3i ^2i,3i 



(66) 



(67) 



with u — cos(z9) and v — sin('i9) e^^ . The minimisation of 
the groundstate energy with 



\HF) ^ a\,a\.al,a\.\ 



(62) 



shows that — for any value of U and the angle d is 
obtained from 



tan4(i?) 



2t 



tan^id) -1 = 0. 



(63) 



The occupation numbers are given by 



1 



(64) 



and shown in Fig. 1151 which illustrates the spontaneous 
symmetry breaking for any value of U. For f7 — > 00 we 
have a perfect antiferromagnet. 

We can now perform a SCRPA calculation in the sym- 
metry broken basis. The RPA operators are given by 



QL = 



^l(T,3<T ^3(T,lcr 



a A— a 



(65) 



In Figs. llBlandlTTIwe give the results. The most striking 
feature is that s-RPA and SCRPA are very close and 
that the error with respect to the exact solution does not 
become greater than 25% for any value of U . Though 
the improvement of SCRPA over s-RPA is very small 
in each channel, at the end in the groundstate energy 
this sums to a more substantial correction in the right 
direction for the ground state energy. This is shown in 
Fig. El a function of atan{Y). We see that HP, s- 
RPA and SCRPA become exact for J7 = and t/ — > 00. 
In between SCRPA deviates e.g. by 8% from the exact 
result at J7 ~ 6 [atan[^) ~ 1.4) whereas this deviation 
is 20% for s-RPA. 

Concluding this section on the 4 -sites case at half 
filling we can say that in the symmetry unbroken basis 
SCRPA is unable to account for some low lying excita- 
tions and therefore fails to reproduce the groundstate en- 
ergy as well. In the symmetry broken basis SCRPA gives 
very little correction over s-RPA. However, the maxi- 
mum error is not greater than 25% for all values of U for 
the excited states and the groundstate energy in SCRPA 
whereas this is 30% for standard RPA. This may be an in- 
teresting result in view of the importance of the so-called 
'plaquettes' (see e.g. ref [23|) in high Tc superconduc- 
tivity. Nevertheless, even though one plaquette (4 -sites) 
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atan(U/t) 

FIG. 16: Energies of excited states with standard RPA, 
SCRPA, and exact solution as a function of U for four sites 
with spin projection = in the symmetry broken basis. 
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FIG. 17: Groundstate energies in HF, standard RPA, 
SCRPA, and exact solution as a function of atan(^) for four 
sites with spin projection = in the symmetry broken 
basis. 



may reasonably be described, the present approach can 
not account for the situation of many plaquettes in inter- 
action which is the real situation in 2D. For the future 
it is therefore very interesting to develope an extension 
of the present SCRPA which not only gives the exact 
solution for the 2 -sites case but equally for the 4 -sites 
case. Such a generalisation is possible in including into 
the RPA operator in addition to the Fermion pair op- 
erators also quadruples of Fermion operators. This is a 
general principle and it has already been demonstrated 
to hold true in the case of the simpler Lipkin model 24| . 
One could call such an extension a second SCRPA in 



analogy to the well known standard second RPA which 
involves in addition to the p/i configurations also 2p — 1h 
ones. In the case of many plaquettes this second SCRPA 
would then constitute a self consistent mean field theory 
for plaquettes. 



V. DISCUSSION, CONCLUSIONS AND 
OUTLOOK 

In this work a many body approach which has essen- 
tially been developed in the nuclear physics context in 
recent years f£| has been applied to the Hubbard model 
for finite number of sites. The theory is an extension 
of standard RPA, called Self Consistent RPA (SCRPA), 
which aims to correct its well known deficiencies as the 
quasiboson approximation with its ensuing violation of 
the Pauli principle and its perturbation theoretical as- 
pect. Of course the appealing features of RPA, as for 
instance fuUfillment of sum rules, restoration of broken 
symmetries, Goldstone theorem, numerical practicabil- 
ity and physical transparency should be kept as much 
as possible. That this is indeed the case with SCRPA 
has in the past been demonstrated with applications to 
several non trivial models (lOj as for instance the many 
level pairing (Richardson) model ^3 and the 3-level Lip- 
kin model |ll|. SCRPA can be derived by minimising an 
energy weighted sum rule and it is therefore a non per- 
turbative variational approach though it is in general not 
of the Raleigh Ritz type. The resulting equations are a 
non linear version of the RPA type which can be inter- 
preted as the mean field equations of interacting quantum 
fluctuations. Though the SCRPA equations are of the 
Schrodinger type, their non linearity non the less makes 
their numerical solution quite demanding. We therefore 
thought it indicated to begin with applications to the 
Hubbard model restricting them to low dimensional cases 
given by a finite number of sites where exact diagonali- 
sation can easily be obtained. We then logically started 
out considering the two sites case (with periodic bound- 
ary conditions), increasing the number of sites by steps 
of two, i.e iV = 2,4,6, .. .. To our satisfaction SCRPA 
solves the 2 -sites problem exactly for any value of U . 
This, as a matter of fact, did not come entirely as a sur- 
prise, since the same happened already with the pairing 
problem for two fermions 10] and indeed it can be shown 
that SCRPA solves a general two body problem exactly 
[l7l |. It is nontheless worth pointing out that other re- 
spectable many body theories fail in the two particle case, 
apart from the low [/-limit. 

In the four sites problem at half filling SCRPA failed. 
This, as in all 4n (n = 1,2,3,.. .) cases, presents the par- 
ticular problem that the system is unstable with respect 
to the formation of staggered magnetisation for any finite 
value of U and this prevented the SCRPA solution to ex- 
ist in the plane wave basis for particular values of the 
momentum transfer \q\. At the end of the paper we indi- 
cated that extending the present RPA ansatz of pft. pairs 
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to include quadruples of fermion operators can solve not 
only the two electrons but also the four electrons case ex- 
actly. This is particularly interesting in view of the fact 
that the 4 -sites case (plaquette) may be very important 
for the explanation of high Tc superconductivity, in con- 
sidering the many plaquette configurations in 2D [23l |. 
In this work we jumped directly to the six sites problem 
which, as all 2 + An cases, causes no particular difficulties 
in SCRPA, even in the symmetry unbroken basis of plane 
waves. Of course, in the case of 6 -sites, SCRPA is not 
exact any more. However, it is shown that the results 
are still excellent for all quantities considered : excited 
states, groundstate, and occupation numbers. Contrary 
to the two sites case, the SCRPA solutions in the plane 
wave basis cannot be obtained for all values of U. Some- 
where after the point where, as a function of U, the first 
mean field instability shows up, the SCRPA also starts to 
deteriorate and in fact docs not converge any longer. Of- 
ten the mean field critical value of U is by passed by 20% 
up to 50% in SCRPA, still staying excellent. However, 
to go into the strong [/-limit we have to introduce the 
above mentioned quadruple fermion operators or perform 
a SCRPA calculation in the symmetry broken basis |3 . 
Such investigations shall be left for the future. We also 
gave arguments why we think that, going to the > 6 
cases, the precision we found for A^ = 6 will not deteri- 
orate. We therfore think that our formalism will allow 
to find precise results for system sizes where an exact di- 
agonalisation becomes prohibitive. Problems in 2D with 
closed shell configurations probably also can and shall be 
considered with the present formalism. Also, as shown 
in ^3 1 the extension to finite temperatures is possible. 

We also should mention that in this work we neglected 
the so-called scattering terms of the form a|, a^' or a\ a^i , 
that is fermion ph operators where either both indices are 
above or both below the Fermi level. In standard RPA 
those configurations automatically decouple from the ph 
and hp spaces. However, in SCRPA with its rounded 
occupation numbers, there is formally no reason not to 
include them. As a matter of fact, as shown in earlier 
work 0, ^1 , to assure the fuUfiUment of the f- sum rule 
and the restoration of broken symmetries, these scatter- 
ing terms must be taken into account. In the present 
case, as well as in earlier studies, the scattering terms 
seem to be almost linearly dependent with the ordinary 
ph and hp configurations. This fact induced difficulties 
with the iteration procedure, since they correspond to 
very small eigenvalues of the norm matrix. Though we 
do not exclude the possibility that this difficulty could be 
mastered with a more refined numerical algorithm, we fi- 
nally refrained from persuing this effort, since we could 
show that the influence of the scattering terms on the 
results is only on the level of a fraction of percent and 
also the f- sum rule is only violated on this order. 

In short we showed that SCRPA, as in previous models, 
performs excellently in the symmetry unbroken regime of 
the Hubbard model. However, the high [/-limit and the 
An sites cases need further developements. 
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APPENDIX A: PARTICLE-HOLE 
CORRELATION FUNCTIONS 



We give the commutations rules which will be useful 
in the calculation of the correlations functions in the ph 
channel. 



l-M,, 



1 - M, 



-23;,rE('^rQt,+3^rQ.O 



,(A1) 



Thus, the following average values can be calculated 
(commuting the Q's to the right) 



ixrxr-y^yn 



^ '(I'-w)) 



ij 



((1-Af,)(l-M,)) 



-2E 



X"^ X"^ X"^ X"" 



yi^3 yi^2 yfi y^o 

-^z -^z -^l -^l 



(A3) 
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Finally, one can express the correlation function accord- 
ing to the amplitudes RPA, {Mi) and of (MiMj) as 



{M^)) 



^ xpxpx:^-xp - yTyTy^y? 



(A4) 



APPENDIX B: DENSITY-DENSITY 
CORRELATION FUNCTIONS 

Given that this RPA formalism preserves the number 
of particles per spin -a (owing to the fact that the trans- 
formation HF does not break the symmetry of spin), one 
has 



p h 



(Bl) 



and the average value {N„) = N^, = ^ what gives us 



(B2) 



On the other hand, one also has 

p h 

p' h' 

with the average value {N^N^') — N^ + N^i, which gives 
us 

p h p' h' 

^<y' ((X ^ho- - X 

p h 

+N,{{Y,np... -Y^nh'a')) (B4) 

p' h' 

Thus for our case, there is the relation 

((X^^^pT -X"''T)(X"p'i -X"'''i)) = 

p hp' h' 

KY.^npa)-Y.{nh.))=Q (B5) 



pa 



ha 



[1] A.L. Fetter, J.D. Walecka, Quantum Theory of Many- 
Particle Systems, McGraw-Hill, Inc. 1971. 
[2] G.D. Mahan, Many-Particle Physics, Plenum Press, 
New York, 1981. 

[3] F. Furche, Phys. Rev. B 64, (2001) 195120 

[4] G. F. Bertsch, C. Guet, K. Hagino, 'physics/0306058' 
[5] A. K. Kerman and C. Y. Lin, Ann. of Phys. (N.Y.) 241 
(1995) 185; 

A. K. Kerman and C. Y. Lin, Ann. of Phys. (N.Y.) 269 
(1998) 55. 

[6] P. Ring, P. Schuck, The Nuclear Many-Body Problem, 
Springer, Berlin 1980 

[7] G. Baym, L. P. Kadanoff, Phys. Rev 124 (1961) 28; 
G. Baym, Phys. Rev 127 (1962) 1391; 
G. Baym, Phys. Lett. 1 (1962) 242. 

[8] J. P. Blaizot, G. Ripka, Quantum Theory of Finite Sys- 
tems, The MIT Press, Cambridge, 1986. 

[9] P. Schuck, S. Ethofer, Nucl. Phys. A 212 (1973) 269; 
J. Dukelsky and P. Schuck, Nucl. Phys. A 512 (1990) 
466; 

J. Dukelsky, G. Ropke, P. Schuck, Nucl. Phys. A 628 
(1998) 17. 

[10] P. Kriiger, P. Schuck, Europhy Lett. 27 (1994) 395; 

J. G. Hirsch, A. Mariano, J. Dukelsky, P. Schuck, Ann. 



Phys. (NY) 296 (2002) 187; 

A. Storozhenko, P. Schuck, J. Dukelsky, G. Ropke and 

A. Vdovin, Ann. Phys. (NY) 307 (2003) 308-334. 

[11] D. Delion, P. Schuck, J. Dukelsky, submitted in Phy. Rev. 

C (nucl-th/0405002). 
[12] A. Rabhi, R. Bennaceur, G. Chanfray, P. Schuck, Phys. 

Rev. C 66 (2002) 064315. 
[13] D. S. Schafer and P. Schuck, Phys. Rev. B 59, (1999) 

1712-1733. 

[14] F. Aryasetiawan, T. Miyake and K. Terakura, Phys. Rev. 

Lett. 88 (2002) 166401. 
[15] G. Seibold, F. Becca, J. Lorenzana, Phys. Rev. B 67, 

085108 (2003). 

[16] Y.M. Vilk, L. Chen, A.-M.S. Tremblay Phys. Rev. B 49, 
13 267 (1994); ibid. Physica C 235-240, 2235 (1994); 
Y.M. Vilk et A.-M.S. Tremblay, J. Phys. I France 7 
(1997) 1309-1368; 

S. Allen and A.-M.S. Tremblay Phys. Rev. B 64, 075115 
(2001); 

B. Kyung, J. S. Landry, A.-M.S. Tremblay, Phys. Rev. 
B 68, 174502 (2003). 

[17] D. Delion, P. Schuck, to be published. 

[18] H.M. Sommermann. Ann. Phys. 151 (1983), p. 163; 

D. Vautherin and N. Vinh Man, Nucl. Phys. A 422 



17 



(1984) 140; 

N. Vinh Mau and D. Vauthcrin, Nucl. Phys. A 445 

(1985) 245. 

[19] M. Grasso and F. Catara, Phys. Rev. C 63, 014317 
(2000). 

[20] D. Janssen and P. Schuck, Z. Physik A 339 (1991) 34-50. 
[21] F. Catara, G. Piccitto, M. Sambataro and N. Van Giai, 



Phys. Rev. B 54 (1996-11) 17536. 
[22] K. Tauaka and F. Marsigho, Phys. Rev. B 60, 3508 
(1999). 

[23] E. Altman, A. Auerbach, Phys. Rev. B 65, 104508 
(2002). 

[24] A. Storozhenko, P. Schuck, to be pubUshed. 



